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Abstract 

The optimized expansion is used to formulate a systematic approx- 
imation scheme to the probabihty distribution of a stochastic system. 
The first order approximation for the one-dimensional system driven 
by noise in an anharmonic potential is shown to agree well with the 
exact solution of the Fokker-Planck equation. Even for a bistable sys- 
tem the whole period of evolution to equilibrium is correctly described 
at various noise intensities. 



1 Introduction 

The Fokker-Planck (FP) equation is widely used to describe non-equilibrium 
systems in physics, chemistry and biology The stochastic approach con- 
sists in representing the most relevant degrees of freedom of the system by 
the variable x driven by noise and deterministic interaction potential U (x, t). 
The time development of the probability distribution W{x^ t) is given by a 



partial differential equation 

r\ r\ r\'2 

g-W{x,t) = LFpW(x,t) := — [U'(x,t)W(x,t) + D^W(x,t) (1) 

where the diffusion coefficient, D, represents a noise intensity and U'{x,t), 
denoting the derivative of the interaction potential with respect to x, is called 
the drift coefficient. The Green's function of the FP equation, P[x,t\x' ,t') 
which fulfils the initial condition P{x,t\x' ,t') = d{x,x'), is called the transi- 
tion probability (conditional probability), since it describes the evolution of 
the probability density from time t' to t: 

W{x,t) = J P{x,t\x',t')W{x',t')dx'. (2) 

For a time independent potential, U{x), the separation ansatz 

W{x,t) = ^{x)e"'' (3) 



L,p*(x) = ^^i±i3i±i + D^ = -K*(x). (4) 



reduces the time dependent FP equation( |T]) to the stationary eigenfunction 
equation 

aU'(x)$(x) j^92#(x) 
9x 5x^ 

The lowest eigenvalue of a FP operator is identically zero, k = 0, and the cor- 
responding eigenfunction ^o{x) can be found exactly, yielding the stationary 
probability distribution of the form 

W,t{x) = %{x) = Ne-^ (5) 

with the normalization constant N = {J^^e~^~dx)~^. For an arbitrary 
potential U{x) the higher eigenfunctions and the non- stationary probability 
distribution cannot be found exactly, there is thus a need to develop approx- 
imation methods. To this end it is convenient to transform the FP equation 
into the Schrodinger equation The transformation 

^(x) = e^$(x) (6) 
brings the FP operator to the Hermitian form 

U(x) -U(x) dp , , 

C = e—Lppe^ = D— - V{x) (7) 
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where 

r(£)]!_£!M. (8) 

The operator hC has the same form as the negative Hamiltonian operator 
for the quantum mechanical particle of the mass M = ^ in the potential 
hV{x). The transformed FP equation 

= (^D^ - V{x)^ ^(x) = A^(a;) (9) 

is called a pseudo-Schrodinger equation, because by (§) the wave function, 
evolves in imaginary time ts = —iht. The transition probability, P{x, t\x' , t'), 
being the Green's function of the original FP equation (|1|), is related to the 
imaginary time evolution amplitude, K{x, t; x', t'), of the pseudo-Schrodinger 
equation by 

P{x, t\x', t') = exp K{x, t; x' , f). (10) 

The evolution amplitude can be represented as a path integral 

K{x,t\x',t')= Dxe-it'^^""^' (11) 

over all functions which begin at x{t') = x' and end at x{t) = x, where the 
pseudo-Schrodinger Lagrangian of a particle is given by 

L[x] = ^ + V{x). (12) 
We shall study a stochastic system in an anharmonic potential, 

2 

U{x) = ^ + Xx\ (13) 
In this case the pseudo-Schrodinger potential (H) takes the form: 
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For vanishing A the problem reduces to the exactly solvable Ornstein-Uhlen- 

2 

beck process in a quadratic interaction potential, U{x) = In this case 
the pseudo-Schrodinger potential is also quadratic, V{x) = — ^ + ^x'^, and 
the path integral for the evolution amplitude (|TTp can be performed exactly 
yielding 



X exp 
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AD sinh(7(t - 1')) 



27rDfl - e-27(*-*')) 

{x^ + x'^) cosh(7(t -t'))- 2xx'] \. (15) 



By d^) this leads to an exact expression for the transition probability of the 
Ornstein-Uhlenbeck process 

-•■^ ■ ' ' ' V2,rO(l-e-2*-''))^l. 2D(i _ e-2T(M')) J • 

In the presence of anharmonicity (A 7^ 0) the evolution amplitude, K{x^ t; x', t') 
cannot be obtained exactly, but various approximation methods are devel- 
oped. The perturbative calculations of P{x,t\x' ,t') in powers of A are pos- 
sible if 7 > 0. In this case the first order approximation describes well the 
evolution of the system approaching the stationary distribution in the long 
time limit, only the normalization worsens with increasing A. However, for 
7 < the method becomes inapplicable, since the transition probability is 
non-normalizable. In this case the perturbative approximations give a wrong 
description of the time evolution, since the maxima of the transition proba- 
bility escape to x = ±00. 

Few years ago we proposed the optimized expansion (OE) scheme for 
the evolution amplitude p| which has much better convergence properties 
than the perturbative expansion in powers of A. Here we apply this method 
to generate a systematic approximation scheme for the transition amplitude 
P(x, t\x' , t') by (0). We will show that the method can be successfully applied 
for describing the time evolution of a stochastic systems. The efficiency of 
the method will be shown on the example of bistable system driven by noise 
in the double-well potential. Such a system attracts much attention in non- 
linear optics (statistical properties of laser light above threshold), solid-state 
physics and chemistry. The solution of the FP equation which describes the 
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evolution from an intrinsically unstable state to the final stationary state 
is of a special interest, and various approximation methods have been pro- 
posed and compared with the numerical solution in this case. However, in 
all the methods the evolution process is divided into few steps and different 
approximations are used in each time sector [Q. We will show that the first 
order result of the OE provides a simple approximation which agrees well 
with the exact solution of the FP equation in the whole period of evolution 
to equilibrium at various noise intensities. 



2 The optimized expansion 

The optimized expansion (OE) has been formulated as a method to generate 
non-perturbative approximations for the effective action in quantum field 
theory |Q . The method consists in calculating the effective action as a series 
in e, by splitting the Lagrangian into 

L = Lo + e(L-Lo), (17) 

where the unperturbed part contains arbitrary parameters, which are op- 
timized in every order calculation (upon setting e = 1). The method is 
equivalent to a systematic re-summation of the perturbation series and gives 
the Hartree-Fock-Bogolubov approximation in the leading order. A similar 
idea has been applied to formulate systematic approximation methods for 
other physical quantities in a number of works P, P] under different names 
(self-similar perturbation theory, delta expansion, variational perturbation 
theory, optimized perturbation theory...). The approach provides a method 
to systematically improve any self-consistent approximation of the theory. 
The OE has been also applied to the quantum mechanics of the particle in 
the potential V{x) by modifying the classical Lagrangian to the form 

L[x]=L^ + ey,„i=^+^— + efy(x) , (18) 

where the harmonic oscillator of the mass M and an arbitrary frequency 
ijj is chosen as the unperturbed system, and calculating physical quantities 
as a series in powers of e. It has been shown, that for the ground state 
energy of the quantum mechanical anharmonic oscillator a convergent series 



5 



is obtained [|^ as opposed to perturbative series which is asymptotic. 
Applying the method to calculate the imaginary time evolution amplitude 0, 
1^, it is convenient to represent the series for W{x,t; x' ,t') = In K{x,t; x'f) 
by the cumulant expansion. For the particle of the mass ^ we have 

W{x,t;x',t') = W^{x,t;x',t')-e <Vintix) >^ 

+ -^i< VLi^) >u;- < Vintix) >^< Vintix) >^) - ...(19) 

where by (|T3p 

W^{x, x', t, t') = In K^{x, x', t, f) = 
to \ 




uo [{x"^ + x'"^) coshijj(t — t') — 2xx'] 



2 \47rDsinha;(t - t')/ 4D sinh a; (t - t' 

and the expectation values are calculated for the unperturbed Lagrangian 



< ... >. 



= / Dx...e-/*'^"["1^*. (21) 

JxH'Ui;' 



The A^-th order approximant, W^^\x't';x,t), is obtained by truncating the 
series (P^ after the A-th term and setting e = 1, since only in this case does 



the modified action agree with the classical one. The exact result, being a 
sum of an infinite series, would not depend on arbitrary frequency, but any 
finite order truncation shows such a dependence. We choose, therefore, the 
value of the unperturbed frequency, u, to make the given order approximant 
as insensitive as possible to its small variation, by requiring 

The optimization condition determines the value of a; as a function of (3, x and 
x', which changes from order to order, improving the convergence properties 
of the approximation scheme. 

In the case of polynomial potential the expectation values in W'^^^ (p!9|) 
are given by Gaussian functional integrals which can be easily performed, 
yielding an analytic expression for W^^\ The first order result for the evo- 
lution amplitude of the quartic oscillator has been obtained in Ref. [@], and 
shown to give a satisfactory approximation to the particle density in the 
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broad range of the oscillator parameters. Here we have to calculate the evo- 
lution amplitude in the potential 

V{x) = go + g2X^ + g4X^ + gax^, (23) 

because the pseudo-Schrodinger potential corresponding to the bistable FP 
equation ([l^ is sextic. After setting e = 1 the first order results can be 
written as 

W^'\x, t; x', t')) = -tgo + W^{x, t; x' , t')+(^^-g2^ /j^'^^^ + K{r)]dr 

-^4 [\L\T) + 6L\T)K{r) + 3K\T)]dr 
Jt' 

-ge [ [L%T) + 15K{T)L\T)+A5K\T)L{Ty + 15K%T)]dT (24) 
Jt' 

where 

^. , xsmh.Lu(T — t') + XbSmh.u(t — T) , , 2Dsinhcij(r — t') sinhcc;(f — r) 
L[T) = . — and K[t) — 



sinha'(t — t') Ci;sinhci;(t — t') 

and the optimization condition (UT) reduces to 



(^?2-y)^/j^^(r) + Alr)]rfr+^4^/j^^(r) + 6L^(r)ir(r) + 3ir2(r)]dr 
+g6^ I\l\t) + lhK{T)L\T) + AhK\T)L'{T) + lhK\T)]dT = 0. (25) 

OUJ Jt' 



Upon performing the integrals in the above expressions an analytic expres- 
sion for the evolution amplitude K{x^t\x' ,t') = e^C^:.*;^; ,* ) jg obtained. The 
result is optimized by choosing the value of the frequency to fulfil (|^) and 
used to calculate the first order approximation to the transition probability 
P{x,t\x' ,t') according to Eq.||. 

After completing this work we learned that approximate solutions of the 
FP equation in an anharmonic potential have been also studied by directly 
improving the perturbative expansion of the transition probability, using 
a drift coefficient as a variational parameter |^. The variational expres- 
sions substantially differ from ours, because the sextic term of the pseudo- 
Schrodinger potential ([T^) does not contribute to the first order in their 
method. The numerical results are also different, especially in the case of the 
bistable potential. 
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3 Results and conclusions 



The approximation obtained in the first order of the OE is exact at A = 0, 
and it gives very good results when the anharmonic terms are small and the 
pseudo-Schrodinger potential is convex (7 > 0). The shape of the transition 
probability is well described, but the total probability slowly decreases in 
time. A similar problem appears in the first order calculation in the vari- 
ational perturbation method 0, where the total probability increases in 
time. The spoiling of the normalization by optimization is a general feature 
of perturbative schemes with variational parameters [0 and can be cured by 
normalizing the each order result by hand. The normalized probability dis- 
tribution obtained in the first order of the OE shows a good agreement with 
the exact solution of the FP equation in the whole period of evolution. Here 
we show the results in the most demanding case of double-well interaction 
potential (7 < 0), when the pseudo-Schrodinger potential has a multiple- well 

2 

structure. The critical value of the diffusion coefficient (Dcr = 24\)^ distin- 
guishes two cases: the pseudo-Schrodinger potential of double-well shape in 
the case of large noise {D > Dcr), and of triple-well shape in the case of small 
noise. 

We compare the results of our approximation to the transition probabil- 
ity, P{x,t\0,0), with the exact results calculated numerically in Ref. ^ and 
Ref. for the interaction potential U{x) = —\x^ + \x'^, which corresponds 
to 7 = —1 and A = I in our notation. The critical value of the diffusion 
coefficient is Dcr = | this case. The results for the time evolution of the 
transition probability at noise intensity D = 0.1 are shown in Fig.l, and the 
results at D = 0.05 and D = 0.01, in Fig.2 and Fig.3, respectively. The 
pseudo-Schrodinger potential for the considered values of D is also plotted. 

The transition probability for the system being initially in the unstable 
state (P(x,0|0,0) = 6{x — 0)) is presented at different times of the evolu- 
tion: the first value is in the initial time region, the next ones are in the 
intermediary region, and the last one is in the final time region when the 
stationary distribution is already achieved. It is remarkable that a good de- 
scription of the evolution from the unstable configuration to the stable one 
is obtained already in the first order of the OE, even in the difficult case of 
small noise. This is due to the optimization of the variational parameter u 
at a given point x in the considered time t. One has to stress, that in our 
calculation the stationarity condition, which determines the optimal value 
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of uj (|25|), has a solution for all the values of x during the whole period of 
evolution. This is in difference with the variational perturbation method [^], 
where different criteria of the variational parameter fixing are used in the 
initial and in the final regions, and the worst description is obtained in the 
intermediate stage of evolution. In our approximation the initial and inter- 
mediate stage of evolution are described very well for all values of noise, but 
the discrepancies appear in the final stage of evolution and the asymptotic 
distribution differs from the exact one. For D = 0.05 and D = 0.01 the 
approximation is of similar quality as the two or three-stage approximations 
based on fl expansion 0, and gives also a good description for D = 0.1, 
when the approximation discussed by Hu becomes inaccurate as observed in 



Ref. |TT|. The accuracy of our approximation can be improved by higher 
order calculation in a systematic way. 

One has to note that an extension of this approach to higher dimen- 
sional systems is possible. Also the dynamics of a bistable stochastic system 
driven by time- dependent forces can be studied in the OE. The influence of 
a periodic force on the bistable system, which is a topic of current interest 
because of the phenomenon of stochastic resonance, will be discussed in a 
future publication. 
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(a) 




Figure 1: The pseudo-Schrodinger potential ( a) and the time evolution of 
P{x, t\0, 0) (h) for the bistable potential U{x) = + x"^ at D = 0.1 in the 
first order of the OE {dashed line) at t = 0.3, t = 0.9, t = 1.2 and t = 6., 
compared with the exact results taken from Ref. [|lT| (solid line) 
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